Analysing Lyapunov spectra of chaotic dynamical systems 
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It is shown that the asymptotic spectra of finite-time Lyapunov exponents of a variety of fully 
chaotic dynamical systems can be understood in terms of a statistical analysis. Using random 
matrix theory we derive numerical and in particular analytical results which provide insights into the 
overall behaviour of the Lyapunov exponents particularly for strange attractors. The corresponding 
distributions for the unstable periodic orbits are investigated for comparison. 
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o 

Finite-time Lyapunov exponents represent an important tool for the quantitative description of the geometrical and 
dynamical properties of chaotic dynamical systems 1] and apply therefore to a variety of different physical situations. 
Numerous recent works (see, for example, refs. |^-6]) try to describe the spectrum of these exponents in terms of 
products of random matrices. The majority of these investigations focus on the maximum average Lyapunov exponent 
which can be estimated by taking the proper mean of an ensemble of chaotic trajectories for a sufficiently large time 
interval p. Details of the corresponding distribution of Lyapunov exponents (LE) are thereby not important. In order 
^| ' to analyze strange attractors in more detail it is however crucial to know the spectrum of the finite-time LE which 
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provides valuable information on the structures and properties emerging in phase space pLM. For the asymptotic 
(large p) Lyapunov spectra it is known ^ that, apart from a few exceptions ||, a Gaussian approximation fits the 
behaviour around the maximum very well. Very little is however known with respect to the overall behaviour of the 
i— h • distribution. This is in contrast to the fact that the non-Gaussian tails of the distribution have significant influence 
, Q , | on physical processes Q . 

In the present letter we analyze finite-time Lyapunov spectra for low-dimensional discrete dynamical systems with 
t-H 1 fully developed chaos for the case of large time intervals p > 1, i.e. their asymptotic form. Our main interest is 
twofold. First we explore common features of these spectra beyond their behaviour in the vicinity of their maxima and 
compare them to the corresponding spectra of the unstable periodic orbits (UPOs). It turns out that these features 
can be understood in terms of products of random matrices. Second we investigate the origin of spectral properties 
which depend on the dynamical system. System-dependent characteristics can e.g. be due to invariant structures in 
| phase space like UPOs. 

It is an open question how the random-like features of the chaotic dynamics determine the distribution of the finite- 
time LE. One way to investigate this is to replace certain dynamical quantities of the original system by random ones 
and to study the resulting changes and/or common properties in the spectra of LE. In several recent works Jll| , ^2[ such 
random-like modifications have been suggested. Adopting this methodology here, we will use four different ensembles 
of trajectories leading to different distributions for the corresponding LE. Each ensemble Ei consists of N trajectories 
(typically N w 10 6 ) of a given length p (p > 1) and is given by (i) the distribution of the starting points of the finite 
trajectories and by (ii) the rule applied to generate the trajectory itself. We will then be primarily interested in the 
distributions of the finite-time maximum LE A = ^lnA belonging to these ensembles {Ei}. A is the largest of the 
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absolute values of the eigenvalues belonging to the transfer matrix — TTMj(^) with {xk\k = 1, ...,p} being 

fc=i 

a finite-time trajectory of the ensemble Ei. Mi(xk) is the stability matrix belonging to the (chaotic) dynamical law 
x n +i = F(x n ). The index i of Mi indicates the ensemble Ei according to which the points {xk} of the trajectories 
are determined. It is important to note that all ensembles {Ei} (sec below) use the specific functional form of the 
stability matrix belonging to F but involve trajectories with different degrees of randomness. In the following we 
specify the ensembles {Ei}. 

E\ consists of trajectories of length p obtained via iteration of the chaotic dynamical law F. The corresponding 
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initial conditions are distributed according to the invariant density of the map. This gives the so-called finite-time 
Lyapunov exponent distribution (FTLED) (see ref. jjj and refs. therein). The ensemble E2 consists of trajectories 
which are generated by a random variable distributed according to the invariant density of the chaotic map F. 
This yields the bootstrap Lyapunov exponent distribution (BLED) , . Compared to the FTLED the BLED 
corresponds to a dynamics with enhanced random character. The successive points of the bootstrap trajectory are 
completely uncorrelated. The third ensemble £3 uses a uniformly distributed random variable for the generation 
of the trajectories. The range of the uniform distribution is chosen according to the phase space of the dynamical 
system. This case corresponds to a random matrix simulation of the dynamical system however respecting the form 
of the stability matrix M which belongs to the map F. The resulting distribution of the Lyapunov exponents is called 
the random matrix Lyapunov exponent distribution (RMLED). Within the present investigation the ensemble E3 
possesses the highest degree of randomness. The fourth ensemble E4 consists of the UPOs of the dynamical system 
F with period p and correspondingly the distribution of their maximal LE. For fixed p the number of trajectories 
contained in E4 is finite according to the topological entropy of the corresponding phase space. 

Let us now explore the distributions defined above for a variety of low-dimensional fully chaotic systems. Our main 
goal is to analyse and understand the overall behaviour of the FTLED for these systems. We begin with a simple Id 
example: the logistic map. Results on the FTLED for this system can be found in ref. Noteworthy the FTLED 
has a non-Gaussian form with one dominating central cusp. In comparison to this our numerical calculations on both 
the BLED and RMLED show that they are smooth functions with a Gaussian-like maximum but with characteristic 
asymmetric tails. For the maximum Lyapunov exponent distribution of the UPOs the exact result is a ^-function i.e. 
p p {\) = S(X — In 2) independent of the period p. Looking at the FTLED one observes that the tails of the distribution 
are rather similar to the tails of the BLED and the RMLED while the cusp (maximum) is located exactly at A = In 2. 
The BLED (or RMLED) reflects therefore the overall behaviour of the FTLED i.e. describes the envelope of the 
FTLED. The latter possesses an additional central peak at the position of the Lyapunov exponents of the UPO's. 
These features will in the following turn out to be common for a broad class of dynamical systems. 

We focus in the following on two-dimensional systems with a strange attractor. The Henon map for a — 1.4 
and b — 0.3 and the Ikeda map jl4| for a = 0.9 and 6 = 6 are prototypes of such systems. The FTLED of both are 
presented in Fig. 1. The ensemble E\ consists of 10 6 trajectories of length p = 27 for both the Henon and Ikeda 
map. The results are surprisingly similar to the Id case: A characteristic envelope with asymmetric tails dominates 
the distribution while superimposed peaks indicate the presence of additional structures in phase space. As we shall 
see in the following three regions with different functional forms of the envelope can be distinguished in the FTLED 
as well as the RMLED and BLED: A fast asymptotic algebraic decay for large values of A, a Gaussian-like behaviour 
around the maximum and a dominating exponential decay for sufficiently small A. Figure 2 shows the BLED and 
RMLED for these maps for the same length p. The position of the maximum of the RMLED is sensitive with respect 
to the random number intervals chosen, i.e. its location carries the information of the position of the attractor in 
phase space. Based on the above results and observations we are naturally lead to the following conclusion: the basic 
possible features of the smooth envelope of the FTLED (asymmetric structure, asymptotic tail properties) of a chaotic 
dynamical system are of random origin and can be obtained and understood by a corresponding study of random 
matrices (see below). Additional superimposed structures are signatures of e.g. invariant sets in phase space and are 
therefore of exclusively deterministic dynamical origin. 

To elucidate and quantify the above observations we perform in the following an analytical investigation of the RM- 
LED. This will allow us to thoroughly understand the behaviour of the RMLED and consequently the corresponding 
aspects of the FTLED. We begin by introducing a fictitious dynamical system with a stability matrix Mj of strongly 
random character i.e. 

where Ti is a random variable uniformly distributed in [0, R] and i labels the fictitious trajectory of length p. The 

simple form of the matrix Mi allows us to factorize the random variables {ri\i — l,...,p} of the transfer matrix 

p 

M^ = J^J Mk and to reduce the problem of the product of random matrices to that of a product of random 

fe=i 

numbers. The matrix structure is then retained in the constant matrix A which is assumed to be nonsingular. 
The distribution of the maximum LE for trajectories of length p of this system is determined as: 

rR rR pR 1 P P 

p p {\) = / dn / dr 2 ... / dr p S(X - - lnTTri|A max |) T[p(n) (1) 
Jo Jo Jo P 
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where p(z) — <d(z)Q(R — z)^, 6 being the step function and A max is given by A max — j(TrA + 
sign(TrA)-^/ (TrA) 2 — 4detA). Using the substitution ti = ln(rj|A m£ia ;|) and performing the Fourier transform of 
the ^-function involved in (pJ) we obtain: 



P P W = ^ I dk I dt v ... I dtp exp [ (1 + ) (2) 



27r(i?|A maa; |)P 



ln(i?.|A maa; |) ,ln(i?|A ma:E |) 



The integrations over ti in (En can be easily performed leading to a p-th order pole in the complex fc-plane at k = i. 
This pole structure, corresponding to values of k for which the exponent of the second exponential term in (Q) vanishes 
and for which the integration of U leads to singularities, is responsible for the features of the LED described above. 
Complex integration finally yields: 

P P W = --^-y (lnGRlA^I) - A)"" 1 e -P(KB|A™.„|)-A) e(ln(i?|A ma;c |) - A) (3) 

Eq.(^|) demonstrates that the exponential behaviour dominates for sufficiently small values of A. Around the maximum 
at A = ln(i?|A maa; |) — (1 — |) the saddle point approximation (p ^> 1) gives a Gaussian. For values of A close to the 
maximum value X rnax = ln(i?|A moa; |) we arrive at a power law behaviour, i.e. an algebraic decay. 

Although the fictitious dynamical model discussed above captures the main features of the statistical properties 
of the distributions of Lyapunov exponents it is clearly desirable to investigate chaotic dynamical systems for which 
the RMLED can be obtained analytically. To this end let us consider the dynamical system defined by the quadratic 
equations: x n+ i = ay£ + b ; y n = cx n + d. This system possesses for a = c = d = 1 and b = —2.5 a strange attractor 
which contains repeating crosses of decreasing size. For that reason we call it in the following the cross map. The 
average maximum Lyapunov exponent of this attractor is A rj 0.123 while its fractal dimension is dp ks 1.78. The 
RMLED for the cross map can be calculated analytically following the line described above for our fictitious model. 
One pecularity of the cross map is that one has to distinguish between the RMLED obtained through trajectories 
with even and odd length p. The reason is that for odd p both eigenvalues of the transfer matrix have the same 
absolute value, while for even p there are two eigenvalues different with respect to their absolute value and one has 
to select the maximal one. After a tedious calculation using complex contour integration techniques we find for the 
RMLED of the cross map for odd values of p the result: 



ft, (A) 
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exp (-p(A -2\))jr(p)(A- J -\n\^\- 2X)^ 9(A - 3 - In & - 2A) (4) 



Here we have taken the random variables appearing in the stability matrix of the cross map to be uniformly distributed 
in the interval [i?i,i?2]- The parameter A is given as: A = ln|2aci?2|- From eq.(Q) we see that the RMLED of the 
cross map is essentially a product of a single exponential and a sum of power laws. The latter possess all the same 
power p — 1 and differ only with respect to the constants involved. The similarity to the RMLED result of our model 
system in eq.(||) is obvious which confirms the universality of certain features of the RMLED. The RMLED for the 
case of even p is given by the integral: 

Pp {\) = 2p»{\) ( dzpp(z) (5) 

J — oo 

with pp(x) according to eq.([|). The integration in (||) can be performed analytically leading to a lengthy expression 
which will be not given here. The main characteristics of the function p p (X) are again the features stated previously: 
a dominating exponential behaviour for sufficiently small values of A, a Gaussian maximum and a fast algebraic decay 
for A close to its maximum value. We have also studied the FTLED, BLED and the distribution of the Lyapunov 
exponents of the UPOs for the cross map for p — 28. The results are shown in Fig. 3. The envelope of these 
distributions exhibit the same features as discussed above. The additional structures present in the FTLED of the 
cross map are, as can be seen from Fig. 3, due to the presence of the UPOs which provide signatures of deterministic 
dynamical origin. 

Finally let us consider the distributions of the Lyapunov exponents of the UPO's of the Henon (p = 27) and the 
Ikeda (p — 14) map Jl6|,[l7j which are presented in Fig. 4. The signs of the characteristic properties of the envelopes 
discussed above are also visible here. An interesting feature appears for the distribution of the LE of the UPOs of the 
Ikeda map: it is shifted significantly compared to the FTLED. This shift is probably due to the fact that the UPOs 
fail to reproduce the invariant density of the attractor, at least up to the above-considered period. It is well-known 
that the Ikeda attractor needs a description going beyond the linear neighbourhood pa|. 
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Summarizing our results we have demonstrated that the overall behaviour of the finite-time Lyapunov exponent 
distributions of fully chaotic dynamical systems show general characteristics i.e. they can be understood in terms of 
statistical random matrix simulations of the systems. Seemingly this holds also for the distributions of the Lyapunov 
exponents of the unstable periodic orbits embedded into the chaotic phase space. Since Lyapunov spectra are at the 
heart of our understanding of chaotic systems in general our results apply to a variety of different physical systems. 
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Figure Captions 

FIG. 1. Fhe FFLED (see text) of the Henon and Ikeda map for p = 27. 
FIG. 2. The BLED and RMLED (see text) of the Henon (-H) and Ikeda (-1) map for p = 27. 

FIG. 3. The FTLED, BLED and the distribution of the Lyapunov exponents of the unstable periodic orbits for the cross 
map for p — 28. 

FIG. 4. The distribution of the Lyapunov exponents of the unstable periodic orbits of the Henon (p = 27) and Ikeda map 
(p=14). 
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Fig. 1 
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